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Abstract 



I The characteristics of the chiral phase transition are analyzed within the framework of chiral 

CO . quark models with nonlocal interactions in the mean field approximation. In the chiral limit, we 

develop a semi-analytic framework which allows us to explicitly determine the phase transition 
curve, the position of the critical points, some relevant critical exponents, etc. For the case 
I of finite current quark masses, we show the behavior of various thermodynamical and chiral 

\^ . response functions across the phase transition. 
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g ; PACS numbers: 12.39.Ki,11.30.Rd,12.38.Mh 

^ ■ I Introduction 

The behavior of strongly interacting matter under extreme conditions of temperature and/or density 
^ . has important consequences in nuclear and particle physics as well as in astrophysics and cosmology. 
^ I From the theoretical point of view, even if a significant progress has been made on the development 
of ab initio calculations as lattice QCD ^ |21 EIj these are not yet able to provide a detailed 
knowledge of the full QCD phase diagram, and most theoretical approaches rely in the study of 
low energy effective models. Qualitatively we expect that chiral symmetry, which is broken at very 
low temperatures and densities, will be restored as the temperature and/or density are increased. 
However, the precise characteristics of this phase transition are still not known. For two massless 
flavors most effective approaches to QCD suggest the existence of a tricritical point on the (T, /j,) 
plane, which separates a first order phase transition line found at lower T and larger fj,, and a 
second order transition line where the chiral restoration occurs for higher T and lower For two 
light flavors a similar behavior is expected, although the second order transition line is replaced 
by a more or less sharp crossover and, correspondingly, the tricritical point becomes an end point. 
In any case, for a given effective model that can provide a reasonable successful description of low 
energy strong interactions, it is very important to obtain as much information as possible about 
the characteristics of the chiral restoration transition. In previous works 0| Ej we have begun the 
study of chiral restoration in the context of chiral quark models with nonlocal interactions [S], 



which can be considered as some nonlocal extension of the widely studied Nambu— Jona-Lasinio 
model j3- In fact, nonlocality arises naturally in the context of several successful approaches to 
low-energy quark dynamics as, for example, the instanton liquid model [HI and the Schwinger-Dyson 
resummation techniques [S]. It has been also argued that nonlocal covariant extensions of the NJL 
model have several advantages over the local scheme like e.g. a natural regularization scheme which 
automatically preserves the anomalies JHI) small NLO corrections [llj . etc. In addition, it has 
been argued |12| I13j that a proper choice of the nonlocal regulator and the model parameters can 
lead to some form of quark confinement, in the sense that the effective quark propagator has no 
poles at real energies. 

Several studies |131 1141 115j have shown that these nonlocal chiral quark models provide a satis- 
factory description of the hadron properties at zero temperature and density. The aim of the present 
work is to complement the analysis of Refs. presenting further details about the chiral phase 

transition within these schemes. Indeed, we show that in the chiral limit it is possible to develop 
a semi-analytic framework which allows us to explicitly determine the phase transition curve, the 
position of the critical points, some relevant critical exponents, etc. For the case of finite current 
quark masses, we present the behavior of various thermodynamical and chiral response functions 
across the phase transition. In particular, it is found that thermal and chiral susceptibilities show 
clear peaks which allow to define a phase transition temperature as it is usually done in lattice 
calculations. 

The paper is organized as follows. In Sect. II we provide a short description of the model and its 
treatment in the mean field approximation (MFA). In Sect. Ill we study the phase transition in the 
chiral limit by performing the Landau expansion of the free energy. We also obtain the MFA critical 
exponents. In Sect. IV we present and discuss the behavior of the different thermodynamical and 
chiral response functions for the case of finite current quark masses, and the corresponding phase 
diagrams are shown. In Sect. V we present a summary of our main results and conclusions. Finally, 
some details of the calculations are given in Appendices A and B. 



II Nonlocal chiral quark models 



Let us begin by stating the Euclidean action for the nonlocal chiral quark model in the case of two 
light flavors and SU{2) isospin symmetry. One has ^ 



Se = d X 



G 

1p{x) {-i^ + mct)lpix) - —ja{x)ja{x) 



(1) 



where tp = {u d)^ and rric stands for the u and d current quark mass. The current ja{x) is given by 



3a{x) = / d y d z f{y - x) f{x - z) ipiy) Ta ipiz) , 



(2) 



where Ta = (1, i75f), and the function f{x — y) is a nonlocal regulator. The latter can be translated 
into momentum space, 

d'^p 



f{x — z) 



(2vr)^ 



<x-z)p _ 



(3) 



^For simplicity we neglect here possible diquark channels. See Ref.|lfci| for details on their role in the phase diagram 
of this type of models. 
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In fact, dimensional analysis together with Lorentz invariance implies that r(p) can only be a func- 
tion of p^/A^, where A is a cutoff parameter describing the range of the nonlocality in momentum 
space. Hence we will use for the Fourier transform of the regulator the form r\{p'^) from now on. 

From the Euclidean action in Eq. the partition function for the model at zero T and /i is 
defined as 

VxpVip e"-^^ . (4) 



We perform now a standard bosonization of the theory, introducing the sigma and pion meson 
fields Ma{x) = ((t(x), vf (x)). In this way the partition function can be written as [H] 

VaVTT det A{Ma) exp -—J ^ M^{p) , (5) 

where the operator A reads in momentum space 

A{Ma) = {-^ + m,) (27r)4 J^^) {p - p') + r^{p^) M,(p - p') r^{p'^) F, . (6) 

In what follows we will work within the mean field approximation, in which the meson fields 
are expanded around their translational invariant vacuum expectation values 

a{x) = a+ 6a{x) (7) 
TTi{x) = 5Tri{x) (8) 

and the fluctuations 6a{x) and S7ri{x) are neglected (vacuum expectation values of the pion fields 
vanish owing to parity conservation). Within this approximation the determinant in Q is formally 
given by 

tr log [-^ + mc + a- r\{p^) ] , (9) 

where tr stands for the trace over the Dirac, flavor and color indices, and V^^'^ is the four-dimensional 
volume of the path integral. 



Ill Phase transition in the chiral Umit 

In Refs. 1^ |5j we have analyzed the chiral restoration within nonlocal chiral quark models for some 
definite regulators. In particular, we have shown that in the chiral limit this phase transition is a 
second order one for low values of the chemical potential ^. In this Section we will rederive this 
result following the so-called classical approach proposed by Landau, in which one expands the free 
energy in powers of the order parameter (in this case the quark condensate {qq)) in the vicinity 
of the critical temperature. This shows the equivalence of the chiral restoration in nonlocal chiral 
quark models with the corresponding phase transitions taking place in other physical systems, such 
as ferromagnets or superfluids. As stated, we will work within the mean field approximation, which 
in this context means to approximate the path integral in Eq. @ by its maximum, reached at some 
saddle point. 

In our case, the partition function in the grand canonical ensemble for finite temperature T and 
chemical potential ^ can be obtained from Eqs. © and ©. In these expressions, the integrals over 
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four-momentum space have to be replaced by Matsubara sums according to 



where ujn are the Matsubara frequencies corresponding to fermionic modes, uJn = (2n+l)7rT. In the 
same way the volume V^^^ is replaced by V/T, V being the three-dimensional volume in coordinate 
space. As in Refs. we are assuming here that the quark interactions only depend on the 

temperature and chemical potential through the arguments of the regulators. The grand canonical 
thermodynamical potential per unit volume is thus given by ^ 

T 

UMF{T,fi,mc) = -— log Zup{T,fi,mc) 



g -4iV,^log[/ + s2(/)] , (11) 



where S(p^) = mc + ar'^{p'^) stands for the quark selfenergy, and the mean field value a{T,n,mc 
is obtained from the condition 

rif,j 

. (12) 



da 

In fact, (Jmf turns out to be divergent. The regularization procedure used here amounts to define 
^m'P'Xt, mc) = uJmf{T, h, rUc) - ujfree{T, /i, nic) + UJ^I'^^}{T, n, rUc) + ujq , (13) 

where wj'^^ (T, fi, TTic) is the regularized expression for the thermodynamical potential of a free 

fermion gas, and ujq is a constant fixed by the condition Lv[fp^^ = Oatr = ^ = (see Appendix A 
for details). 

In the analogy with a ferromagnetic system, the chiral condensate can be identified with the 
magnetization per unit volume, {qq) < — > —M/V, whereas the current quark mass rric plays the 
role of the external magnetic field H. One has 

where the ^ factor results from the fact that this relation holds for each quark flavor separately, 
i.e. q = u,d, while they have a common mass rric- In the chiral limit, the existence of a second 
order phase transition for a fixed value of implies that the condensate (qq) goes to zero when the 
temperature T approaches from below a given critical value Tc{fi), above which one has (qq) = 
and the chiral symmetry is restored. Thus, for T Tc{fM) and to leading order in mc, one can 
perform the Landau expansion (see Appendix A) 

JrF'\T,fi,mc) = ujo + Jj';£{T,fi,mc = 0) + A{T, fi) {qq)^ + C{T,fi){qq)^ + 

+ 2mc{qq) + 0{{qq)Mqqf m^, ml) , (15) 
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where the coefficients A and C are given by 



A{T,fi) 
C{T,fi) 



±-N,S2i{T,f,) 



with 



l28N^Sf,{T,fi) 32 N^SUT,fi) 



— -iVe52l(r,/i) 



Smn 



,2n 



(16) 



(17) 



The regularized thermodynamical potential for a massless fermion gas — second term on the r.h.s. 
of Eq. ((T3|l — can be calculated by evaluating the integral in Eq. HA.2|) in the massless case. One 
obtains 



CO 



?;/i(r,/.,o) = -ii£ 



30 ^ 27r2 ^ 



1 



(18) 



In the limit rric = 0, it can be seen |18j that for C > the system undergoes a second order 
phase transition at a critical temperature Tc{fi) obeying A{Tc{fj.) , fi) = 0. This implies 

1 



(19) 



which defines a second order transition curve in the (T, ^u) plane. As it is described in Appendix B, 
for sufficiently small (but relevant) values of T and /i the Matsubara sum implicit in 52i(Tc(/i), /x) 
can be analytically worked out. This leads to a simple relation between Tc and /i, namely 

/9o-?^^, (20) 



A 



A^,GA2 ' 



where (3q = A"^ J dp p rj^{p'^) is a dimensionless quantity which depends only on the shape of the 
regulator. Note that this relation generalizes that obtained in Ref. for the NJL model. We also 
point out that for a given value of ^u, and temperatures which are close to the critical value T(.{^) 
obtained from Eq. H19|) . one has 

A(T,^l) = \t, (21) 



where 



t^{T-Up:))/T,{^l) , 



(22) 



48 NcSf,{Tc{fi),fi) 

We note in passing that once the Landau expansion has been established, it is a usual textbook 
exercise jl8| to derive the critical exponents ruling the behavior of the specific heat, the order 
parameter (qq) and the chiral susceptibility near the second order critical points. 



Here the specific heat and the chiral susceptibility are defined as 

'd'^uj{T,fi,'mc 



Xv,T,ii 



1 f d'^u!{T, fi, rric 
2 



^ ^ /'ds{T,fi,mc 



dT 

f d{qq){T,fi,mc) 
\ drric 



(23) 

(24) 
(25) 
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where s(T, n, nic) is the entropy density. As expected, one obtains the mean field critical exponents^ 



a = 



/3 



1 

2 ' 



7 



1 



(26) 



a result that was under discussion in the somewhat related Dyson-Schwinger models of QCD |19j . 
Eq. ()26() implies that the chiral susceptibility diverges at T = Tc. In the case of the specific heat, 
there is no divergence but a finite discontinuity at rric = 0, t = 0, 



^v,ti,mc\t=0- 



-v,ti,mc \t=0+ 



+ 



NcT^ 



36 S, 



42 



(27) 



where c 

(reg) 



v,fj,,mc\t-- 



^0+ is the specific heat at constant due to the free fermion gas pressure pfr 



-u; 



free 



{T,fj.,0), evaluated at T = Tc. From Eq. ifTH]) one gets 



-V,fi,mc\t>0 



(28) 



As stated, the transition remains second order as long as C {Tc{fJ,) , fj,) > 0. This is expected 
to be the case for low values of the chemical potential. However, if /x is increased, one could 
reach a point at which the coefficient C vanishes. Beyond this point, called "tricritical" , the 
system undergoes a first order phase transition (the order parameter (qq) is discontinuous at the 
corresponding critical temperature). In this way, the tricritical point is obtained from the conditions 
A{T, fi) = C{T, /i) = 0, or equivalently 



5'42(rc(/u),/x) = . 



(29) 



Again, for sufficiently small (but relevant) values of T and fi the Matsubara sum in can be 
analytically worked out (see Appendix B), leading to 



where 



Pi 
P2 



vr 



/3i (T/AY - + if^/ry ]+P2- fif^/T) - log(r/A) 



(30) 



8 A 



2 dr^ip'^) 



dp 
P L 



dp^ 



p=0 



+ 2(7-1) -log 2vr , 



(31) 



and f{x), which satisfies /(O) = 0, is given by Eq. (|B.10|) of Appendix B. It is interesting to notice 
the similarity between our expressions in Eqs. (|2()j) and (|29ll3()j) and those obtained in Ref. [201 
within a very different theoretical approach. A brief discussion on the subject is also included in 
Appendix B. 

Let us consider for definiteness a nonlocal model in which the regulator is a Gaussian function, 

rA(p2) = exp( V/2A2) . (32) 



^At the tricritical point, the mean field critical exponents are a — ^, /3=j,^— 1 and 5 = 5. 
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In this case one has Pi = 4, while the integral in P2 vanishes. In addition, it is easily seen that 
Sa2{T, 0) is a positive function of T/A, which implies that in the chiral limit the model leads to a 
second order chiral phase transition for vanishing chemical potential. The transition remains second 
order when ^ is increased up to the tricritical point, and the corresponding transition line in the 
{T, fj,) plane can be immediately obtained from Eq. (|2flj) . with /3o = 1/4. The critical temperature 
at ;U = is thus given by 

V3 , A 4vr2 



(33) 



We illustrate these features by considering a particular parameter set, namely A = 760 MeV and 
G = 30 GeV~^. These are the parameters corresponding to Set II in Ref. [5], leading to Tc{0) = 102 
MeV. The second order phase transition is clearly shown in Fig. 1(a), where the solid line represents 
the chiral condensate for /i = as a function of T. Now using Eq. ()29() one can find the position 
of the tricritical point, which in this case is found to be located at {T,fi) =(72 MeV, 133 MeV). In 
Fig. 1(a), dashed and dashed-dotted lines show the behavior of the chiral condensate for fi = 100 
MeV and /i = 200 MeV, below and above the tricritical point respectively. Both the values for 
Tc{0) and the position of the tricritical point are found to be in very good agreement with those 
numerically obtained in Ref. [S]. Finally, the described effect on the specific heat is shown in Fig. 
1(b), where we plot Cy^^ as a function of the temperature for /_f = 0, 100 MeV and 200 MeV. 
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Figure 1: (a) Chiral condensate vs. temperature, and (b) specific heat Cy,^ vs. temperature. Solid, 
dashed and dashed-dotted lines correspond to /i = 0, /i = 100 MeV and /i = 200 MeV respectively. 

The result in Eq. (|33j) is useful to point out a generic feature of this type of models in the 
chiral limit. In fact, it is not difficult to show that, in such limit, for zero T and /i one obtains 
expressions for the pion decay constant fi^^ and the chiral condensate such that fi^^ /A 

and (qq)^^") /A depend only on the combination GA^. In order to get a ratio f^^/{qq)^^^ — 0.4, as 
required by phenomenology, these expressions lead to Gh? ~ 16. If this value is now replaced in 
Eq. one gets rc(0)/A ~ 0.12, which is somewhat large, but still lies within the expected range 
of validity. However, if one imposes the critical temperature T'c(O) to be approximately equal to 
170 MeV, as suggested by lattice QCD calculations, one gets A ~ 1.4 GeV, which, together with 
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Gh? ~ 16, enhances f^x^ and (qq)^'^^ up to roughly 40% above the corresponding phenomenological 
values. Thus, it can be generically said that, if in the chiral limit — and within the mean field 
approximation — one wants to satisfy the phenomenological constraints on the values of /tt and 
(qq), the nonlocal models under consideration lead to a relatively low critical temperature Tc{0). 
Of course, this might be modified by finite quark mass effects and/or beyond MFA corrections. 



IV Phase transition and response functions for finite quark mass 

In this Section we concentrate in the analysis of the phase transition for finite current quark masses 
in the isospin limit. For this purpose we will consider different thermodynamical response functions. 
In particular, it will be seen that, for reasonable values of the current quark mass rric, the chiral 
susceptibilities show clear peaks that can be used to define the transition curves in the crossover 
region. 

The response functions are basically given by the second derivatives of the Helmholtz free energy 
F(T, V, N). The latter is related with the thermodynamical potential by 

F{V,T,N)=[niV,T,i,) + 2N , (34) 

where N = {d^/d^)v,T (as in the case of the quark condensate, the ^ factor comes from 
defining N = Nu = N^). We will consider as before the limit of a large system in which ^{V, T, fi) = 
Vu;{T,fi), F{V,T,N) = V f{T, N/V), thus instead of average particle number N our quantities 
will be given in terms of the average particle density p = N/V. One can define three independent 
thermodynamical response functions, namely the specific heat at fixed volume and particle number 
Cv,p,m, the isothermal compressibility HT.p.mc and the coefficient of thermal expansion Op^p^m^. These 
are given by 



T fd^FiT,V,N,mc)\ ^fds{T,p,mc 



-V,p,mc 



^ O^FiTVN,m.)\ ^fdpiTp,m.)\ ^^^^ 



^T,p,m, V dV^ JT,N,m, V dp y ^^^^ 

= HT,p,m, 7^ (37) 



There are also other thermodynamical response functions that can be defined, such as the specific 
heat at constant pressure Cp^p^m^ and the adiabatic compressibility Hs^p^m^- However, these can be 
written in terms of the three quantities in Eqs. (|35I37|) . 

Notice that we have explicitly included the dependence of the free energy on the current quark 
mass rUc- In fact, one can still define some extra response functions associated with the chiral 
transition by differentiating the free energy with respect to rric- In the analogy with magnetic 
systems, these would be the specific heats at constant magnetization and constant applied magnetic 
field, Cmi Ch, the isothermal and adiabatic susceptibilities xt, Xs, and the coefficient of thermal 
magnetization, ajj- As before, these quantities are related to each other, leaving only two new 
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independent response functions. Here we choose to consider those analogous to the isothermal 
susceptibility and the coefficient of thermal magnetization, thus we define 

1 ( d^F{T,V,N,m,) \ _ ( d{qq){T,p,mc) \ 



1 (d^F{T,V,N,mc)\ fd{qq)iT,p,m, 



2V V dm^dT Jy^^ V dT , ^^^^ 



(39) 



Let us show the numerical results for these quantities for the case of a nonlocal Gaussian 
regulator. We have chosen the parameter set A = 760 MeV, G = 30 GeV~^ and nic = 7.7 MeV, 
which corresponds to Set II in the notation of Ref. It can be seen that these parameters lead 
to the empirical values of the pion mass and decay constant, and provide reasonable results for the 
chiral quark condensate and the quark selfenergy S(0) at zero T and /i. In Fig. 2 we show the curves 
corresponding to the specific heat Cv^p^mc ™d the chiral response functions Xv,T,p and av,p,mc- Iii ttie 
left panels, these functions are plotted versus the temperature for three representative values of the 
density (the quoted values are referred to nuclear matter density, po — 1-3 x 10^ MeV^). Notice that 
for low nonzero densities the system is homogeneous only for temperatures which exceed a critical 
value. Below this limit, as we will see, there is a region where phases with broken and restored chiral 
symmetry coexist. In the right panels, the same response functions are plotted as functions of the 
density, fixing the temperature at 50 MeV and 100 MeV. Once again, for low temperatures there 
is a mixed phase region, and as a consequence the functions are not well defined at intermediate 
densities. We have also analyzed the response functions corresponding to the parameter set given 
by G = 50 GeV-2, mc = 10.5 MeV and A = 627 MeV, or Set I in the notation of Ref. 0. Set I and 
II might be interpreted as confining and nonconfining respectively, where confinement is understood 
in the sense that the pole structure of the quark propagator does not allow quarks to materialize 
on-shell in Minkowski space |12[ I13j . The curves in the case of Set I do not differ qualitatively from 
those shown in Fig. 2, therefore they have not been included here. 

Now let us pay special attention to the responses of the chiral condensate (order parameter of the 
phase transition) at fixed temperature and chiral quark mass (or "magnetization"). These are given 
by Eqs. (|38|) and (|39j) and their behavior as functions of the temperature is shown in the second 
and third rows of Fig. 2 (left panels). As it is well known, susceptibilities are particularly useful 
to analyze the phase transition features. In the present model, as shown in previous works 
for low temperatures the system undergoes a first order chiral phase transition, which turns into 
a smooth crossover for temperatures exceeding a given "end point". In this crossover region, the 
transition is characterized by the presence of respective peaks in the mentioned susceptibilities, the 
height and sharpness of these peaks giving a measure of the crossover steepness. The same can be 
observed if one looks at the chiral and thermal susceptibilities for fixed /i, which are given by the 
second derivatives of the thermodynamical potential, and are in fact the natural quantities to deal 
with in the grand canonical ensemble. The definition of the chiral susceptibility Xv,T,p has been 
already given in Eq. (|25|) . while the thermal susceptibility at constant fi is defined as 

1 f dMT,P,m,) \ f d{qq){T,fi,m,) \ 
av,p,m, = - o ( ) = - ( 7^ ) ■ (40) 

/ p,mc 



2 V dnicdT /, V dT 



The curves showing the behavior of the susceptibilities Xv,T,p and av,p,m as functions of the 
temperature and the chemical potential are shown in Fig. 3 (the chosen parameter set is the same 
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as in Fig. 2). For completeness we also include the curves for the specific heat at constant /i, 
previously introduced in Sect. Ill (where the case of varying rric was discussed). In a natural way, 
the peaks in the curves for Xv,T,fi and av,^,mc can be used to define the position at which the 
chiral transition occurs. Thus, one can extend the phase space diagram to include the crossover 
transition curves in addition to the first order ones. This is represented in Fig. 4, where crossover 
curves obtained from the peaks in Xv,T,fi and av,fi,mc have been represented by dashed and dotted 
lines respectively. Although chiral and thermal susceptibilities lead to slightly different points, it 
is seen that the transition region is well defined. Full lines correspond to the first order phase 
transition, while the fat dots indicate the end points. Left and right panels show the results for 
parameter Sets I and II respectively, while upper (lower) panels show the phase space transition 
curves in the T — fi [T — p) plane. Notice that in the T — p phase diagrams there is a region below 
the first order transition lines where both phases are allowed. The latter can be interpreted [22 
as a zone in which droplets containing light quarks of mass nic coexist with a gas of constituent, 
massive quarks. The dotted lines inside this zone are the spinodals, i.e. the boundaries of the region 
in which the energetically unfavored solutions can exist as metastable states. 

Our results for the phase transition curves in the T — p plane are qualitatively similar to those 
obtained from lattice QCD calculations |21 E], although both the critical temperature at ;U = 
and the end point temperature turn out to be relatively low in our case (see discussion at the 
end of Sect. III). For the sake of comparison, it is also interesting to analyze the curvature of the 
phase boundary at p = 0. This is an appropriate quantity to be studied in lattice QCD, where the 
main problem in the analysis of the T — p phase diagram is the inclusion of a finite real chemical 
potential. Recent lattice calculations P yield Tc{(PTc/dp'^)\^=o = —0.14 it 0.06, while the result 
obtained within the standard NJL model (up to some finite current quark mass corrections) gives 
a value of about —0.40 j22j . In our model, the curvature can be calculated numerically from the 
results plotted in Fig. 4, leading to values of about —0.26 and —0.30 for Sets I and II, respectively. 
In the chiral limit the corresponding calculation can be carried out from the analytical expression 
in Eq. leading to Tc{cPTc/ dp'^)\^=o = -S/vr^ ~ -0.304. 

Finally, let us quote the numerical results for the behavior of the energy density, the entropy 
density and the pressure as functions of the temperature. The corresponding curves are shown in 
Fig. 5, where we plot the scaled quantities e/T^, | s/T^ and 'ip/T'^ versus the relative temperature 
T/Tc, for Sets I and II and different values of the chemical potential. The arrows indicate the 
value corresponding to a free fermion system in the large T limit, given by 77rVl0 —see Eq. (ITKl) . 
Notice that for Set I there is a range of temperatures below Tc for which all three quantities are 
negative. This might be taken as a further indication that for this parameter set there is a sort 
of "confinement" : somewhat below hadronic degrees of freedom ought to be included, and the 
system cannot be simply treated as a quark gas. This represents a qualitative difference between 
Sets I and II. It is interesting to see that, in spite of this fact, the behavior at the transition region 
is relatively similar in both cases. 

V Summary and outlook 

In this work we have presented further details of the chiral phase transition within chiral quark 
models including nonlocal interactions in the mean field approximation. In the chiral limit, our 
analysis allows to obtain semi-analytical expressions for the transition curve and the location of 
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the critical points. For the case of finite current quark masses, we have studied the behavior of 
various thermodynamical and chiral response functions across the phase transition. In the crossover 
region the thermal and chiral susceptibilities display clear peaks which allow to define a transition 
temperature, in the same way as it is usually done in lattice calculations. 

The resulting phase diagrams are qualitatively similar to those obtained within other frame- 
works, as e.g. the standard NJL model, the Bag model, and lattice QCD. In the nonlocal schemes 
studied here, however, the transition temperature at = is found to be somewhat lower than the 
expected value Tc(0) ~ 170 MeV once the model parameters are fitted so as to reproduce both the 
empirical values of and /tt and the phenomenological value of the light quark condensates at 
zero temperature and density. In fact, this result may be improved by considering alternative ways 
of including nonlocality, as e.g. regulator schemes inspired in one-gluon exchange processes |23j . or 
through the inclusion of beyond MFA contributions. Work in this direction is currently in progress. 
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Appendix A: Expansion of the thermodynamical potential in pow- 
ers of {qq) 

The grand canonical thermodynamical potential defined by Eq. Hll|) turns out to be divergent. Here 
it has been regularized by subtracting the expression corresponding to a system of free fermions of 
mass rric, and adding it in a regularized form. We have thus 



(T, /i, rric 



a 

2G 



4NrM ^^ log 



+ ml 



where 



■ANrT 



(2^ 



log (^1 + e-(^-^)/^) + log (l + e 



(A.l) 



(A.2) 



with E = Y^p2 + m§ , while ujq is a constant fixed by imposing that the thermodynamical potential 
vanishes at zero T and ^. We point out that this regularization prescription relies on the well- 
behaved shape of the regulator in S(p^). Alternative schemes can be applied in other cases, see 
e.g. El. 

Let us assume that in the chiral limit, rric = 0, the order parameter (qq) vanishes at a given 
temperature T = Tc- We are interested in the description of the situation in the vicinity of this 
temperature, hence it is useful to carry out a double expansion of oj^mf^'^ i^i powers of {qq) and rric- 



For rric = 0, the expansion of Uf,fp 
derivatives 



{Teg) ■ 



in powers of {qq) can be performed by taking into account the 
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^ gV/^ (d{qq)y' d^qq) 
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(qq)=0 
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d{qq)" 



(qq)=0 



(A.3) 



Since {qq) = implies o" = 0, it is easy to see that the above partial derivatives are given by 



d{qq)_ 
da 



da"^ 



a=0 



a=0 



-4Nc SniT,fi) 

^ - 8N,S2i{T,fi) 



where the functions Smn(T,fj,) are defined by 
On the other hand, for nic 7^ one has 

duJMF 



d^jqq) 
da^ 

da^ 



a=0 
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,2n 
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dnir 



2{qq) 
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mc=0 



Taking into account this relation, together with Eqs. ()A.4() and HA.3|) . one easily arrives at the 
result shown in Eq. H15|) . 



Appendix B: Evaluation of Matsubara sums 

In this Appendix we will show how to work out the Matsubara sums Smi and Sm2 1 which are defined 
by Eq. (ICT . These sums appear in the Landau expansion of the thermodynamical potential near 
the critical temperature, see Eq. (fTK]) . 

In order to carry out the calculations, we take into account the analysis in Ref. where 
Cauchy's theorem was used to convert the Matsubara sums into an integral plus a sum over pole 
residues. In the case of a function -F(p^) having only simple poles and no cuts in the complex plane 
one obtains ^ 

j—^ Re 2J 7fc Res [Tiz); Zk] [n+{zk + /u) + n_(zfc + 1^1)] , (B.l) 
4 >o 

where the function T is defined as T{z) = F[{—iz — i^Y + IpI^Ij = Rk + Hk are the residues of 
this function in the complex plane z, and the coefficient 7^ is defined as 7^ = 1/2 for 1^ = and 
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7A: = 1 otherwise. We have also introduced here (complex) occupation number functions n±{z), 
defined by 

' 1 + exp [(2 ^ /i)/r] ^ ^ 

In the case of the sum 52i(T, //), the corresponding function J^2i{z) is given by 



^21(2) 



-(z + ^)2 + |p12 



(B.3) 



which has only two simple poles located at = —ji it with residues =p(2|p|)~^. Therefore, 
Eq. (|B.1|) can in principle be applied. There is, however, an subtle point to be taken into account. 
The derivation of Eq. HB.1|) assumes that when Re z — > 00, which is in fact true for all 

the situations considered in Ref. jSj. However, depending on the specific form of the regulator, this 
condition might be not satisfied by the function !F2i{z). If this is the case, Eq. ()B.H) can still be 
applied provided that T,fi<^A (see below). 
Assuming that Eq. (|B.1|) holds, one gets 



S2i{T, fi) 
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p2 
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dp P riip"^) 
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The last integral in this expression can be worked out analytically. One has 

^2 







dp p [ n+ (p) + n_ (p) 



+ Li, 



6 



+ 



(B.5) 



which together with Eq. (|B.4|) leads to the result quoted in Eq. ()2U() . As we have mentioned, this 
relations are only valid for sufficiently low values of T, fj, compared with the cutoff scale A. In the 
case of the Gaussian regulator we have checked that for T < A/6 and /x < A/4 Eq. (|2Ujl is verified 
with an accuracy of less that 1%. In the case of the Lorentzian regulator the region of validity is 
somewhat smaller, but at the same time the relevant cutoff parameter is larger than in the Gaussian 
case As a conclusion, we find that in all cases considered here Eq. H2(J|) can be taken to be valid 
with very good approximation for the values of T and /x of physical interest. 

A similar analysis can be carried out for the sum 5*42 (T,/i). Here the situation is somewhat 
more complicated, since one has to deal with double poles and Eq. (jB.lj) is no longer valid. One 
has instead 
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Res 



1 + e^/T 
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1 + e-^/^ 



, (B.6) 



where F{p^) = r\{p'^)/p'^, ^{z) = F{—{z + fi'^) + Evaluating the residues, this leads to 



-00 



- [n+{p) +n^{p)' 
p 



(p2 + g2^2 

4m r'j^iO) p[n+{p) +n-{p)] - ^ [n+{p) + n-{p) - n^ip) - niip)] 



(B.7) 
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where r^(0) = dr\{p'^) /dp^\p2^Q . The first two terms lead to divergent integrals that can be worked 
out with the help of a definite regulator, e.g. the Gaussian one, writing 



For the Gaussian regulator, the sum of the divergent integrals in (|B.7|) . properly regularized, is 
given by 

g-4{p2+g2)/A2\^ 1 ^ 

-[n+{p) + n-{p)]) = 



where 




3o {p^ + q^Y j P 
1 + 2 _ log(^V^) - log (T/A) - /(^/T) , (B.9) 



f{x) = 2 sinh2(x/2) r^-^ t.nHy/2) 

Jq y cosh X + cosh y 

Finally, the integral of the last term in Eq. ()B.7|) can be explicitly performed, giving simply 



/•oo 

/ dp[n+ip) + n-{p)-nl{p)-nl{p)] = T . (B.ll) 
Jo 

The previous result in Eq. (|B.5() . together with Eqs. (|B.7() to (|B.11() . lead to the expression quoted 
in Eq. (|3U() . Again this relation is strictly valid for T,fi <^ A. For the Gaussian regulator we have 
checked that it still holds up to 1% in the region T < A/9, /i < A/5. As before, this is also the 
case for the Lorentzian regulator for the values of T and n relevant for our analysis. 

It is interesting to point out the similarity between our equations determining the second order 
transition line and the tricritical point, and those obtained in Ref. [SOI- In that article, the authors 
address the study of the T — ^ phase diagram using a different approach, in which they propose 
a large flavor number expansion and a resummed renormalization scheme. It can be seen that the 
results in Ref. [20] — Eqs. (13) and (14) — have the same form as our expressions if one takes the 
limit A = 0, which in the context of Ref. [20] means to neglect the contributions from the meson 
sector (this should be analogous to the MFA considered here). Moreover, the validity of these 
results is also limited to relatively low values of T and in view of the presence of a Landau pole 
arising from the renormalization procedure. 
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Figure 2: Some response functions for fixed p, corresponding to parameter Set II. In the left panels 
curves are given as functions of T, with pjp^ fixed at three representative values, namely 0.25 
(dotted), 0.75 (dashed) and 2.0 (solid). Right panels show the curves as functions of p/po for T 
fixed at 50 MeV (sohd) and 100 MeV (dashed). 
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Figure 3: Some response functions for fixed /i, corresponding to parameter Set II. In the left panels 
curves are given as functions of T, with fj, fixed at three representative values, namely 100 MeV 
(dotted), 250 MeV (dashed) and 300 MeV (solid). Right panels show the curves as functions of /x 
for T fixed at 50 MeV (sohd) and 100 MeV (dashed). 
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Figure 4: Phase diagram showing the chiral transition curves in the T — fi plane (upper panels) 
and the T — p plane (lower panels), for parameter Sets I (left) and II (right). Solid lines correspond 
to first order phase transition curves, while dashed (dashed-dotted) lines show the crossover curves 
obtained from the peaks in the chiral (thermal) susceptibilities. The dotted lines in the lower panels 
represent the spinodals, and the fat dots indicate in each case the position of the end point. 
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Figure 5: Numerical results for the scaled quantities ejT^ (solid), |s/T^ (dashed) and 'ip/T^ 
(dashed-dotted). Curves are given as functions of T/Tc, for fixed values of (jl that are indicated 
in each case together with the respective critical temperatures. The arrows show the asymptotic 
values, which correspond to a free fermion system. Left and right panels show the results obtained 
for Sets I and II, respectively. 
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